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Abstract 

Measurements of the IV characteristics of site-diluted Josephson-j unction arrays have revealed 
intriguing effects of percolative disorder on the phase transition and the vortex dynamics in a two- 
dimensional XY system. Different from other types of phase transitions, the Kosterlitz-Thouless 
(KT) transition was eliminated with the introduction of percolative disorder far below the percola- 
tion threshold. Even after the KT order had been removed, the system remained superconducting 
at low temperatures by establishing a different type of order. Near the percolation threshold, evi- 
dence was found that, as a consequence of the underlying fractal structure, the critical dynamics 
of the phase degrees of freedom persisted down to zero temperature. 
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Percolation has been used as an idealized simple model for understanding a variety of 
irregular systems. The effects of such a geometry on phase transitions have been of continu- 
ous interest in condensed matter physics, the Kosterlitz-Thouless (KT) phase transition [l| 
is an essential factor for understanding various systems in two dimensions, such as super- 
conductors, easy-plane ferromagnets, crystals, and liquid crystals, which can be represented 
by the XY model. In addition, the KT transition is singular with regard to excitations of 
topological defects, which drive the transition from a quasi-long-range ordered phase to a 
disordered phase. Nevertheless, the relationship between percolative disorder and the two- 
dimensional (2D) phase transition has not been sufficiently explored and remains unsettled 

as 

. Also, of interest are the effects of the underlying fractal geometry on the dynamics 
of the phase degrees of freedom of XY systems near the percolation threshold. In this 
paper, we report an experimental study of the effects of percolative disorder on both the 
2D phase transition and the critical dynamics in proximity-coupled Josephson-j unction ar- 
rays (JJA's), which provide a near-ideal realization of the XY model. Measurements of the 
current- volt age (/V) characteristics demonstrate that the nature of the phase transition in 
the 2D XY system is completely altered with progressive addition of percolative disorder. 
Near the percolation threshold, experimental evidence suggests that the underlying fractal 
structure results in dynamical criticality being maintained down to zero temperature. 

Experiments were performed on a series of site-diluted square arrays of Nb/Cu/Nb 
Josephson junctions. Nb crosses were periodically disposed on a Cu film with a lattice 
constant of 13.7 jttm, a junction width of 4 /im, and a junction separation of 1.4 fim to form 
a 400x600 square lattice. Nb islands were randomly removed from the lattice to introduce 
percolative disorder. Arrays with the concentration of filled sites p = 1, 0.94, 0.86, 0.8, 0.75, 
0.7, 0.65, and 0.6 were fabricated. For the samples, the same random-number seed was used 
to select the sites to be removed. The percolation-threshold concentration p c for the series of 
samples, below which no percolating cluster of connected junctions was formed, was 0.5906 

a. 

The IV characteristics were measured by employing the phase-sensitive volt age- signal- 
detection method with a lock-in voltmeter and a square- wave current at 23 Hz. The single- 
junction critical current i c and the single-j unction coupling energy J (= hi c /2e) were deter- 
mined by using the de Gennes formula e 5] for proximity-coupled junctions in the dirty limit 
to extrapolate the i c vs T data at low temperatures and by using the numerically obtained / 
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vs dV/dl data for site-diluted arrays in Ref. [6( to compensate the systematic error arising 
from the presence of diluted sites. The i c vs T data at low temperatures were obtained 
from the I vs dV/dl measurements. The magnetic field or the frustration /, defined as the 
number of flux quanta per unit cell, was adjusted by using the R vs f curve of the sample 
exhibiting distinct resistance minima at integral /'s. The temperature during the measure- 



ments was contro 
are given in Ref. 



led to have fluctuations ^ 1 mK. Additional details of the measurements 

a. 

Figure ^ presents some of the results of the IV characteristics measurements on the 
samples of p = 1, 0.86, 0.7, and 0.65, with the ambient magnetic fields expelled from the 
sample space by the /i-metal shield and solenoid. The measurements were carried out at 
17-19 different temperatures (at intervals shorter than shown in the figure). The IV curves 
were obtained by averaging 15-240 measurements for each current. It is visible that the IV 
characteristics change completely with decreasing p. Unlike the p = 1 (or no site disorder) 
sample experiencing a KT transition [8|, the samples with p = 0.7 and 0.65 exhibit low- 
temperature IV characteristics with an activated character. The variation of the voltage 
signal with temperature both above and below the transition is significantly slower for p = 
0.7 and 0.65 than for p = 1 jjj. For p = 0.86, the IV characteristics are more like those 
of a KT system at T > T c , but they exhibit an activated character at T < T c . The 

IV characteristics for p = 0.7 and 0.65 are never similar to those of a KT system jioj . 
The low-temperature curves for p = 0.7 and 0.65 can be fitted to an exponential form 

V ~ Jexp[— (It/I)^] with fx = 0.9-1.1. The low-temperature curves of an exponential 
form indicate that at strong site disorder, the arrays have genuine superconductivity with 
long-range phase coherence. The d(logV) / d(\ogI) vs I/I c plots for p = 0.7 and 0.65 in 
Fig. Q showing the proposed criterion for a non-KT-type superconducting transition to 
be satisfied, confirm a finite-temperature non-KT-type superconducting transition for the 
strongly disordered samples. The evolution of the TV characteristics with p demonstrates 
that percolative disorder far below the percolation threshold alters the nature of the 2D 
phase transition in an unfrustrated XY system from a KT-type transition to a non-KT- 
type transition. 

Figure El shows the IV characteristics of the same samples with frustration / = 2/5. The 
IV plots for / = 2/5 demonstrate that at strong disorder or at lower concentrations of filled 
sites, the IV characteristics are affected little by the imposition of frustration, excluding 
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the /-dependent phase-transition temperature. It is surprising that at strong disorder, the 
IV characteristics of the unfrustrated system become similar to those of frustrated systems, 
which are known to experience a melting transition of a vortex solid driven by domain 
wall excitations Q]. The development of the unfrustrated system into a system similar 
to frustrated systems with the addition of percolative disorder is also revealed in the p 
dependence of the reduced phase-transition temperature T c (= T C /{J /ks)), shown in Fig. 
EJ The T c 's in the figure are the temperatures at which / ~ V 3 for systems exhibiting KT- 
like IV characteristics or at which the IV curve becomes straight in a log /-log V plot for 
systems exhibiting non-KT-like IV characteristics Q]. The T c 's determined from the IV 
curves agree, within the error bars, with those independently determined from the resistance 
measurements. The T c vs p plot shows that the p dependences of T c for / = and 2/5, 
despite the total difference at weak disorder, become qualitatively similar at strong disorder. 
The rapid drop of T c to zero in the vicinity of p c is a common feature of systems with a gap 
in the excitation spectrum The close resemblance of the IV characteristics and the p 
dependence of T c at strong disorder for two systems, one with and one without frustration, 
appears to suggest that the non-KT-type low-temperature order in an unfrustrated XY 
system at strong disorder is possibly that of a superconducting vortex solid with long-range 
phase coherence, as in frustrated systems. 

Although a vortex solid in an unfrustrated array sounds peculiar, the possibility of its 
existence is notable because an unfrustrated array with the KT transition suppressed always 
contains sufficient unbound vortices to form a vortex solid. A recent numerical study [15j 
has suggested that in the presence of a random potential above a critical disorder strength, a 
breaking of ergodicity due to a large energy barrier against vortex motion may allow unbound 
vortices to form a vortex solid at finite temperatures. In the equivalent interacting-vortex 
(or Coulomb-gas) representation of the XY model, random site dilutions produce a random 
pinning potential acting on vortices. Our speculation is that the random pinning potential 
eliminates the KT transition far below the percolation threshold and constrains the unbound 
vortices to form a superconducting solid at low temperatures jlfil ]. 

Another important feature of the data presented in this paper is that the IV charac- 
teristics of the unfrustrated system completely change again near p c (= 0.5906). As Fig. |U 
reveals, the low-temperature IV curves shift from exponential-type curves to power-law-type 
curves when p changes from 0.65 to 0.6. For p = 0.6, an approximate power-law character 
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is maintained even at such low temperatures as T = 0.006, far below T c (= 0.16) |l7|. The 
power-law character never weakens, even when frustration is turned on. In addition, the 
slope of the isotherm at T c in the log /-log V plot is only 1.9 for p = 0.6, less than it is for 
p ^ 0.65. Since the slope minus one is the dynamic critical exponent [l^J, the lower slope 
for p = 0.6 implies that the critical relaxation in a percolating system is faster than it is in 
less-disordered systems. This is contrary to the general expectation that more diluted sites 
in a JJA usually result in slower relaxation of the system. Even if the power-law IV rela- 
tion is characteristic of a KT transition, the / independence of the power-law behavior and 
the dynamic exponent being much less than 2 indicate that the possibility of a re-emerging 
KT transition near p c is low. The motion of the vortices created at diluted sites by an 
injected current does not explain the observed power-law behavior either. If the vortices 
created at the defects were responsible for the power-law behavior of the p = 0.6 system, 
as expected from numerical studies [3| of diluted arrays at zero temperature, a similar be- 
havior should be visible in systems with comparable defect concentrations. The exponential 
low-temperature TV curves for p = 0.65, however, confirm that the contribution of vortices 
created at defects must be insignificantly small. 

Near the percolation threshold, a self-similar fractal structure with geometrical inhomo- 
geneities and a divergent percolation correlation length characterize percolating systems. 
Therefore, even at zero temperature where the thermal correlation length, i.e., the spatial 
extent of the thermal fluctuations vanishes, the spatial configuration of the phase angles 
of the superconducting islands in the percolating array has a self-similar critical structure, 
as at critical temperature. For non-KT-type superconducting transitions, a power-law IV 
relation is a salient feature of dynamical criticality The power-law IV curves far below 
the transition for p = 0.6 may then indicate that temperature-independent criticality due to 
the fractal geometry applies to the dynamics of the phase angles, as well as to their spatial 
configuration. This means that the relaxation time near the percolation threshold becomes 
divergent at all temperatures below T c . It has been suggested [19] from measurements of the 
frequency-dependent complex conductance of a site-diluted 2D JJA that the fractal struc- 
ture results in a crossover, at a critical frequency, of the dynamics of the phase angles from 
a low-frequency Euclidean (homogeneous) region dominated by vortices to a high-frequency 
fractal (inhomogeneous) region in which the relevant excitations are localized fractons. The 
critical frequency was found from the ac measurements to be ^ 1 kHz for arrays comparable 
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to the present sample. However, the results of the IV characteristics measurements on the 
present sample appear to indicate that the low-frequency vortex dynamics is not Euclidean 
either, but is in another inhomogeneous region in which the fractal geometry leads to critical 
dynamics. 

To conclude, the IV characteristics measurements on site-diluted JJA's have revealed 
intriguing effects of percolative disorder on the phase transition and the vortex dynamics in 
a 2D XY system. Unlike other types of phase transitions, the KT transition was eliminated 
by introduction of percolative disorder far below the percolation threshold. When the KT 
order was removed, the system unexpectedly remained superconducting by establishing a 
different type of order with long-range phase coherence until the percolation threshold was 
attained. Near the percolation threshold, as a consequence of the underlying self-similar 
geometrical structure, evidence was found for the critical dynamics of the phase degrees of 
freedom persisting down to zero temperature. 
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FIG. 1: Development of the IV characteristics with the filled-site concentration p for / = 0. The / 
vs V isotherms differ by a temperature interval of 0.1 K. The panels on the right show the slopes of 
the isotherms as functions of I/I c - The dashed lines are drawn to show where the phase transition 
occurs. For p = 1, the dashed line represents the IV curve at the temperature where V ~ I 3 . 
The dashed lines for p = 0.7 and 0.65 are drawn at the temperatures where the IV curves are 
estimated to become straight. For p = 0.86, V ~ I 3 at the transition temperature where the IV 
curve becomes straight. See the text for discussions. 

FIG. 2: I vs V and I/I c vs <i(log V)/d(log I) plots for the samples with p =0.7 and 0.65 at / = 
2/5. 

FIG. 3: p dependence of the reduced phase transition temperature T c (= T c / '(J /&b)) for / = 
and 2/5. The arrow points to the percolation threshold p = 0.5906. The inset shows a portion of 
a site-diluted array used in the experiment. 

FIG. 4: IV characteristics of the sample with p = 0.6 for / = and 2/5. The I vs V isotherms 
differ by a temperature interval of 0.1 K. 
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